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Abstract: Hydrothermal condition is mismatched in arid and semi-arid regions, particularly in Central Asia 
(including Kazakhstan, Kyrgyzstan, ‘Tajikistan, Uzbekistan, and Turkmenistan), resulting many 
environmental limitations. In this study, we projected hydrothermal condition in Central Asia based on 
bias-corrected multi-model ensembles (MMEs) from the Coupled Model Intercomparison Project Phase 6 
(CMIP6) under four Shared Socioeconomic Pathway and Representative Concentration Pathway (SSP-RCP) 
scenarios (SSP126 (SSP1-RCP2.6), SSP245 (SSP2-RCP4.5), SSP460 (SSP4-RCP6.0), and SSP585 
(SSP5-RCP8.5)) during 2015—2100. The bias correction and spatial disaggregation, water-thermal product 
index, and sensitivity analysis were used in this study. The results showed that the hydrothermal condition is 
mismatched in the central and southern deserts, whereas the region of Pamir Mountains and Tianshan 
Mountains as well as the northern plains of Kazakhstan showed a matched hydrothermal condition. 
Compared with the historical period, the matched degree of hydrothermal condition improves during 
2046-2075, but degenerates during 2015—2044 and 2076—2100. The change of hydrothermal condition is 
sensitive to precipitation in the northern regions and the maximum temperatures in the southern regions. 
The result suggests that the optimal scenario in Central Asia is SSP126 scenario, while SSP585 scenario 
brings further hydrothermal contradictions. This study provides scientific information for the development 
and sustainable utilization of hydrothermal resources in arid and semi-arid regions under climate change. 
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1 Introduction 


Hydrothermal condition is a comprehensive climate indicator. It is affected by multiple climate 
variables, including precipitation, temperature, and evapotranspiration, and plays a significant role 
in evaluating regional climate change (Weltzin et al., 2003; McCain and Colwell, 2011; Wu et al., 
2011; Konapala et al., 2020). Hydrothermal condition is characterized by contradiction in arid and 
semi-arid regions, particularly in Central Asia (Wang et al., 2010; Li and Ma, 2018; Luo et al., 
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2019; Yu et al., 2020). In the last century, the pronounced hydrothermal contradiction, a mismatch 
between water and thermal resources, was triggered by phased warming and highly variable 
precipitation, which were not synchronized in intensity and magnitude (McCain and Colwell, 
2011; Lacombe et al., 2012; Luo et al., 2019; Konapala et al., 2020), bringing negative impact on 
agricultural yield and water utilization structure and leading to desertification (Zhan et al., 2018; Li 
et al., 2020; Schierhorn et al., 2020). Increasingly noticeable hydrothermal contradictions will 
occur in the future, and the increasing stress brought by population growth and consequential 
agricultural demands will boost these contradictions (Lioubimtseva and Henebry, 2009; Wu et al., 
2011; Li et al., 2020). Consequently, identifying the consistency between potential changes in 
hydrothermal resources has received extensive attention. However, previous studies have not been 
able to account for the comprehensive connection between hydrothermal resources or to project the 
potential changes in future periods (Wang et al., 2010; Zhan et al., 2018; Yu et al., 2020). 
Therefore, it is necessary to project hydrothermal condition on spatial-temporal scales and to 
minimize the possible negative effects of climate change in Central Asia. 

Many comprehensive indices are used to assess hydrothermal condition, and they are divided 
into three types: a direct hydrothermal index, a comprehensive index based on potential 
evapotranspiration (PET), and a comprehensive index based on temperature and water balance 
(Meng et al., 2004). The first type overlooks a comprehensive understanding of water 
consumption, with examples including Winkler index, Selianinov hydrothermal coefficient, and 
De Martonne aridity index (Martonne, 1926; Seljaninov, 1966; Meng et al., 2004; Cardoso et al., 
2019). PET is considered as the natural consumption in the second type, with the indices based on 
drought (Mondal et al., 2021), humidity (Vicente-Serrano et al., 2010), the Palmer drought 
severity index (Alley, 1984), and so on. 

With the comprehensive examination of water balance and thermal properties, new 
hydrothermal indices have evolved from basic indices. For example, according to the relationship 
among temperature, precipitation, PET, and soil water balance, Ni and Zhang (1997) developed 
water-thermal product index (k index) to assess the matched and mismatched hydrothermal 
condition. Relative to other indices, k index considers more reasonable and comprehensive 
hydrothermal condition and has special ecological significance. It has been validated in the 
research of the connection between climate and large-scale vegetation patterns in China (Ni and 
Zhang, 1997; Li et al., 2018). It can be seen that the application of k index is feasible in 
large-scale regions with spatial differences. Thus, it was chosen as a basic index to assess the 
regional hydrothermal condition. 

Global Circulation Models (GCMs) have made excellent advances stemming from the Coupled 
Model Intercomparison Project (CMIP). Together with their simulations, they are regarded as 
effective tools for climate change evaluation (Eyring et al., 2016). The Coupled Model 
Intercomparison Project Phase 6 (CMIP6) has an upgraded physical algorithm that avoids 
overestimating precipitation and equates reliable data sources for climate projection (Gidden et 
al., 2019; Rivera and Arnould, 2020). The basis scenarios in CMIP6 are designed as a 
combination of Shared Socioeconomic Pathways and Representative Concentration Pathways 
(SSP-RCP) scenarios, and are expected to produce more reliable projections (Gidden et al., 2019). 
Varied SSP-RCP scenarios indicate a gradual increase in the emission level of greenhouse gases 
and radiative forcing. 

This study illustrates hydrothermal condition in typical arid and semi-arid regions in future 
periods. First, bias correction and spatial disaggregation (BCSD) method was used to correct the 
simulated bias of climate variables (e.g., precipitation, temperature, maximum temperature 
(Tmax), minimum temperature (Tmin), and PET) from the selected GCMs. Second, k index was 
adopted to assess hydrothermal condition in Central Asia from 2015 to 2100. Additionally, 
sensitivity analysis was applied to identify the sensitivity of hydrothermal condition to climate 
variables. The results provide fundamental data and theoretical support for future development 
strategies in Central Asia with an integrated index for regional climate analysis. 
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2 Materials and methods 


2.1 Study area 


Central Asia (35°09'35"—55°26'33"N, 46°29'57"—87°20'25"E) is located in the center of Eurasia, 
with a total area of about 3.97x10’ km’, including five countries: Kazakhstan, Kyrgyzstan, 
Tajikistan, Uzbekistan, and Turkmenistan (Fig. 1). Most regions of Central Asia have continental 
characteristics, with excess evapotranspiration as well as extreme changes in daily temperature, 
annual temperature, and seasonal precipitation. The region is divided into the northern plains 
(Kazakhskiy Melkosopochnik, Ustyurt Plateau, and West Siberian Plain), the eastern plateaus 
(Tianshan Mountains, Pamir Mountains, and Altay Mountains), the western lowlands (Turan Plain 
and Caspian Depression), and the southern deserts (Karakum Desert and Kyzylkum Desert). The 
major rivers in Central Asia (Ertix River, Amu Darya River, and Syr Darya River) originate from 
the eastern mountains and flow westward into the Aral Sea. The land use patterns in Central Asia 
are limited by continental climate and fragile ecosystems. Sparse vegetation, grasslands, and 
croplands cover 46.09%, 24.26%, and 20.43% of the total land area, respectively (Jiang et al., 
2022). Central Asia enjoys greater feasibility in spatial analysis of hydrothermal condition and 
contradiction because of its complex and diversified natural landscapes. 


Elevation (m) | 
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“Low: -134 — River —”* International boundary 0 20 km 


Fig. 1 Overview of the study area. Note that this map is based on the standard map (No. GS (2016) 1666) of the 
Map Service System (http://bzdt.ch.mnr.gov.cn/) marked by the Ministry of Natural Resources of the People's 
Republic of China, and the base map has not been modified. 


2.2 Climate data 


Datasets used in this study consisted of historical reference data and GCMs simulations. 
Historical reference data were obtained from the Climate Resource Unit (CRU) TS4.04 
(https://data.ceda.ac.uk/badc/cru/), including PET, precipitation, temperature, Tmax, and Tmin 
during 1975-2014, with 0.5°x0.5° spatial resolution. Existing studies have recognized the validity 
of the dataset for climatic analysis in Central Asia (Mannig et al., 2013; Wu et al., 2013; Harris et 
al., 2014; Guo et al., 2021). GCMs obtained from CMIP6 (https://esgf-node.IInl.gov/search/ 
cmip6/) under four SSP-RCP scenarios (SSP126 (SSP1-RCP2.6), SSP245 (SSP2-RCP4.5), 
SSP460 (SSP4-RCP6.0), and SSP585 (SSP5-RCP8.5)) included monthly downwelling and 
upwelling radiations, precipitation, temperature, Tmax, Tmin, near-surface wind speed (NWS) at 
a height of 10 m, and humidity during the historical period (1975-2014) and the future period 
(2015-2100). In consideration of the availability of all climate variables under four SSP-RCP 
scenarios, we utilized five GCMs, including CanESM5, IPSL-CM6A-LR, MIROC6, 
MRI-ESM2-0, and FGOALS-g3, to assess hydrothermal condition. The bilinear interpolation 
method was used to unify GCMs simulations into the same spatial resolution. In addition, 
considering the uncertainty of GCMs simulations, we also calculated the multi-model ensembles 
(MMEs) in this study (Weiland et al., 2012; Guo et al., 2021). 
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2.3 Data analysis 

2.3.1 Taylor diagram and Taylor skill score 

Taylor diagram and Taylor skill score (TSS) are effective tools to assess the simulation level of 
GCMs (Taylor, 2001; Guo et al., 2021). Taylor diagram was chosen to examine the differences in 
performance of climatic reproducibility between individual GCMs and MMEs in Central Asia. It 
integrates three evaluation metrics: Pearson correlation coefficient (r), standard deviation (SD), and 
mean square error (MSE). On Taylor diagram, the closer the distance between the model and 
reference, the higher the simulation capability. TSS has been widely used to assess the suitability of 
a selected model in a given study area. It can be calculated as follows: 


4x (+r)? 


p 
SDhist | SDovs x a Hr y 
SDo5 SDhyist i 


TSS = 4 (1) 


where TSS is the Taylor skill score; r is the Pearson correlation coefficient; SDhist is the standard 
deviation of historical GCMs simulations; SDops is the standard deviation of historical reference 
data; and rp is the ideal correlation coefficient, which is usually 0.999. The higher the TSS, the 
better the simulation capability. 

2.3.2 Bias correction and spatial disaggregation 


The process of bias correction is a prerequisite step to enhance the reliability of GCMs 
simulations (Christensen et al., 2008; Sun et al., 2019). BCSD is a completed pre-processing 
method for GCMs simulations by integrating bias correction and non-parametric transformation of 
downscaling (Wood et al., 2002, 2004; Christensen et al., 2008). Compared to other bias correction 
methods, such as the Delta method, change factor method, and linear scaling (Teutschbein and 
Seibert, 2012), BCSD corrects the inherent and random bias in the simulated mean and variance. It 
also corrects the intensity distribution by establishing the cumulative distribution function and 
transfer function between the historical reference data and GCMs simulations. BCSD is superior to 
the quantile-mapping method, which is a fundamental requirement in climate analysis 
(Chernozhukov et al., 2017). Therefore, BCSD was chosen as a basis method to correct the 
simulated bias in the historical period (1975-2014) and future period (2015-2100). The MMEs, 
rather than individual GCMs, were applied to handle the uncertainty of GCMs simulations. 

BCSD was used to correct the simulated bias in precipitation, temperature, and PET, which are 
the basic variable for the hydrothermal index. It also corrects the bias in an indirect hydrothermal 
index, such as Tmin or Tmax. Other related variables were transformed into PET based on the 
Penman-Monteith equation and were used to complete the bias correction. The main process of bias 
correction includes the following four steps. 

Firstly, the bilinear interpolation method was used to unify GCMs simulations to the same spatial 
resolution before calculating the MMEs (0.5°x0.5° spatial resolution) (Su et al., 2021; Zeng et al., 
2021). Secondly, the simulated trend was extracted from the raw MMEs with the moving average 
method and then the trend was retained. New simulations with no trends were generated. Thirdly, 
the new simulations in the second step were corrected by the following equations: 


Fic-humg = CDF '(CDF(F; hMME |e F 7 MME); obs Oo Pos) , (2) 


Foe-aume = CDF! (CDF (Fields 0 Mme)» Hobs» 0 obs) > (3) 
where F'bc_hmme is the historical MMEs after bias correction from 1975 to 2014; F»c-tmme is the 
future MMEs after bias correction from 2015 to 2100; CDF is the cumulative distribution function; 
Fume is the historical simulations, which was obtained in the second step; Fmme is the future 
simulations, which was also obtained in the second step; 4nmme and Lobs are the location parameter 
of fitting Gaussian distribution based on the historical simulations and historical reference data, 
respectively; and Owe and o*,,, are the variance of fitting Gaussian distribution based on the 
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historical simulations and historical reference data, respectively. Based on the monthly scale, we 
constructed the transfer function for each spatial grid and climate variable (precipitation, 
temperature, Tmax, Tmin, and PET). 

Lastly, the trend was added to the bias correction results and fitting outliers of cumulative 
distribution function was fixed by extrapolation method. 
2.3.3 Evaluation indices of bias correction 


We chose six indices to evaluate the results of bias correction, and the detail description, equation, 
and optimum value of those indices are shown in Table 1 (Seo and Ok, 2013; Ji et al., 2020). 


Table 1 Evaluation indices of bias correction 


Index Description Equation Optimum value 
i=l 
Bias Bias between the historical SMME isi — ODS 0.0 
MMEs and reference data Bias = 
n 
Relative bias The ratio of bias to historical RB= Bias 0.0 
(RB) reference data ‘obsnist 
Time correlation Time correlation between the DOMME a, Pica Di(0bS ia acs 
(TC) historical MMEs and reference TC= = T 1.0 
data [Bote -MME nis)? x [Bow —obshist )” 
i=l i=l 


m 


Spatial correlation between the DAMME sis, — MMEnist) x D (obs yi.) Obst) 


historical MMEs and reference SC = | = z 1.0 


Spatial correlation 
(SC) 


m m 


data (MME)... ; - MMEnist )” x [Bom j — ObSnist)” 
j=l jal 

al : 
Normalized Normalized standard deviation IMME. MME nisni) 
standard between the historical MMEs and NSTD = -———= 0.8-1.2 
deviation (NSTD) reference data S (obs; —obsnist.i)? 

ist i ashe 
i=l 

Mean absolute Mean absolute error of the >|MME vic = ODS pict, 0.0 
error (MAE) historical MMEs MAE =~ i 


n 


Note: i indicates the year; j indicates the spatial grid; n is the total number of year; m is the total number of spatial grid; MMEhis is the 
historical multi-model ensembles (MMEs); MMEpis is the mean value of historical MMEs; obshix is the historical reference data, and 
ObShis:_ iS the mean value of the historical reference data. 


2.3.4 Water-thermal product index (k index) 

According to the positive correlation and principle of hydrothermal balance, Ni and Zhang (1997) 
developed a regional climate indicator named water-thermal product index (k index) to denote the 
regional hydrothermal condition and contradiction. The equations are as follows: 


k =Tx WD +100, (4) 
WD =P-PET, (5) 
0.408A(R,-G)+R X (E,-E,) 
PET = T+273 , (6) 
A+R(1+0.34u, ) 


where k is the water-thermal product index (°C-mm); T is the temperature (°C); WD is the water 
deficit (mm), indicating the difference between precipitation (P; mm) and PET (mm); R, is net 
surface radiation (MJ/(m?.d)); G is the soil heat flux (MJ/(m2.d)); R is the psychrometric constant; Es 
is the saturated water vapor pressure (kPa); Ea is the actual water vapor pressure (kPa); A is the 
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slope of the saturated vapor pressure curve; and u, is the wind speed (m/s). 

The same sign between water deficit and temperature means the matched hydrothermal condition, 
while the contrary sign between them implies the mismatched condition. The recent studies 
suggested that k index is well applied in the agricultural planning and vegetation classification (Ni 
and Zhang, 1997; Li et al., 2018). 

In this study, we divided the monthly scales into spring (from March to May), summer (from June 
to August), autumn (from September to November), and winter (from December to February) to 
assess the seasonal discrepancies in hydrothermal condition. 

2.3.5 Sensitivity analysis 
The sensitivity coefficient (S) was chosen as an analysis method to reveal the sensitivity of k index 
to the climate variables (precipitation, temperature, Tmin, and Tmax) (Beven, 1979). The equation 


is as follows: 
2 Ak / k 
vo ims o[ ? (7) 


where S, is the sensitivity coefficient of k index to climate variables v; and Ak and Av are the 
variation of k and v, respectively. A positive value of S, indicates the consistent change between k 
index and climate variable, while a negative value of S, implies the inconsistent change between k 
index and climate variable. The higher the absolute value of S,, the greater the sensitivity. 


3 Results 


3.1 Evaluation of simulation performance 


The result showed that the spatial pattern of each climate variable can be captured by the five 
GCMs with distinctions in simulation performance. FGOALS-g3 had the smallest spatial bias in 
precipitation simulation, ranging from —2.66 to 1.07 mm/d, and RB was 0.40. IPSL-CM6A-LR 
achieved better reproducibility in temperature simulation, with RB of —0.04. MRI-ESM2-0 had a 
better performance for PET simulation, with RB of 0.04. Moreover, the simulation performance 
of different GCMs in the same region was different (Figs. 2—4). For example, CanESMS5 
overestimated precipitation in the northern Central Asia and underestimated that in the region of 
Pamir and Tianshan mountains; the overestimated precipitation and temperature were presented in 
IMIROC6 and MRI-ESM2-0. Furthermore, the comparison between the raw MMEs and 
individual GCMs showed that RB of precipitation was lower in MMEs, with the mean value of 
0.39. The mean RB of temperature and PET in MMEs was 0.18 and —0.05, respectively. In 
addition, the selected models had better simulation performance in Central Asia, as evidenced by 


(a) CanESMS (b) IPSL-CM6A-LR (c) MIROC6 (d) MRI-ESM2-0 N 


(e) FGOALS-g3 (f) MMEs 


Legend 
Daeg CaN 


RB 


piet: 1 
F 0 1000km 
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Fig. 2 Spatial distribution of relative bias (RB) of precipitation based on five Global Circulation Models 
(GCMs) (CanESM5 (a), IPSL-CM6A-LR (b), MIROC6 (c), MRI-ESM2-0 (d), and FGOALS-g3 (e)) and 
multi-model ensembles (MMEs, f) in Central Asia during 1975-2014. Note that this map is based on the standard 
map (No. GS (2016) 1666) of the Map Service System (http://bzdt.ch.mnr.gov.cn/) marked by the Ministry of 
Natural Resources of the People's Republic of China, and the base map has not been modified. 
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(a) CanESMS (b) IPSL-CM6A-LR (c) MIROC6 (d) MRI-ESM2-0 N 
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Fig. 3 Spatial distribution of RB of temperature based on five GCMs (CanESMS (a), IPSL-CM6A-LR (b), 
MIROC6 (c), MRI-ESM2-0 (d), and FGOALS-g3 (e)) and MMEs (f) in Central Asia during 1975-2014. Note that 
this map is based on the standard map (No. GS (2016) 1666) of the Map Service System (http://bzdt. 
ch.mnr.gov.cn/) marked by the Ministry of Natural Resources of the People's Republic of China, and the base map 
has not been modified. 


(b) IPSL-CM6A-LR (c) MIROC6 (d) MRI-ESM2-0 N 


Fig. 4 Spatial distribution of RB of potential evapotranspiration (PET) based on five GCMs (CanESMS (a), 
IPSL-CM6A-LR (b), MIROC6 (c), MRI-ESM2-0 (d), and FGOALS-g3 (e)) and MMEs (f) in Central Asia during 
1975-2014. Note that this map is based on the standard map (No. GS (2016) 1666) of the Map Service System 
(http://bzdt.ch.mnr.gov.cn/) marked by the Ministry of Natural Resources of the People's Republic of China, and 
the base map has not been modified. 


Taylor diagram (Fig. 5). In particular, the MMEs had a relatively high ranking and large TSS 
value, which showed that the approach was reliable for the analysis of climate change. 

Overall, the MMEs are superior in the climatological analysis compared with individual GCMs 
because their combination of individual output uncertainty generally minimizes spatial-temporal 
bias of multiple meteorological variables, despite accompanied with dipped edges in single or 
several-element simulations. 


3.2 Evaluation of bias correction 


The statistical evaluation of bias correction is shown in Table 2. Smaller bias and RB in each 
corrected climate variable (precipitation, temperature, PET, Tmax, Tmin) indicated better 
performance of BCSD. According to Table 2, for the five climatic variables, time correlation (TC) 
and spatial correlation (SC) have been significantly increased after correction. The increased 
correlation between historical reference data and simulations in MMEs improved the accuracy of 
spatial-temporal analysis. 

Figure 6 revealed the comparison of interannual relative bias before and after correction. The 
interannual relative bias of precipitation after correction was closer to the optimum value than the 
raw precipitation with decreased spatial anomalies, demonstrating that BCSD improved the 
simulation performance for precipitation and minimized interannual bias for PET and temperature. 
Although a slight increase in corrected MAE was found in temperature, the MMEs corrected by 
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Fig. 5 Taylor diagram of precipitation (a), temperature (b), and PET (c) in Central Asia during 1975-2014. The 
green lines represent the reference of mean square error. 


Table 2 Statistical evaluation of bias correction for climate variables 


index Precipitation Temperature PET Tmax Tmin 

BC AC BC AC BC AC BC AC BC AC 
Bias 0.27 0.01 1.20 0.08 —0.37 0.01 0.44 0.00 1.00 0.04 
RB 0.39 0.01 0.18 0.01 0.15 0.00 0.03 0.00 0.63 0.02 
TC 0.36 0.66 0.96 0.96 0.96 0.97 0.97 0.97 0.97 0.96 
SC 0.62 0.63 0.88 0.96 0.85 0.91 0.91 0.96 0.88 0.96 
NSTD 0.86 0.98 1.09 1.00 1.28 1.01 1.15 1.01 1.10 1.01 
MAE 0.47 0.39 2.97 3.15 0.74 0.41 2.99 2.32 2.87 2.03 


Note: PET, potential evapotranspiration; Tmax, maximum temperature; Tmin, minimum temperature; BC, before correction; 
AC, after correction. 


BCSD were generally close to the references for the whole period and region, and were reliable 
for projecting hydrothermal condition. 


3.3 Projected temporal variation of k index at annual and seasonal scales 


Temporal variation of k index at annual and seasonal scales was projected during the periods of 
2015-2045 (near-term), 2046-2075 (mid-term), and 2076-2100 (late-term) under different 
SSP-RCP scenarios (Fig. 7). Under SSP126 scenario, winter had the best hydrothermal condition 
during 2015-2045, with a maximum k index of —1.57°C-mm, followed by spring and autumn, 
with k index being —21.79°C-mm and —20.29°C-mm, respectively. The severe imbalanced 
hydrothermal condition appeared in summer, with a minimum & index of —140.49°C-mm. The 
declining trend of k index was found at annual and seasonal scales. During 2015-2045, k index 
was lower than that in the reference period, indicating that Central Asia may experience the worst 
hydrothermal contradiction. 


202205.00176v1 


chinaXiv 


ChinaXiv@ ERAT! 


YAO Linlin et al.: Projection of hydrothermal condition in Central Asia under... 


‘ 1.0 
Ae (a) Precipitation (b) Temperature 
0.5 0.5 } 
2 0.0 0.0 
-0.5 -0.5 } 
-1.0 -1.0 
1975 1980 1985 1990 1995 2000 2005 2010 2015 1975 1980 1985 1990 1995 2000 2005 2010 2015 
1.0 
(c) PET (d) Tmax 
0.5 + 
9 0.0 0.0 
-0.5 -0.5 } 
-1.0 -1.0 
1975 1980 1985 1990 1995 2000 2005 2010 2015 1975 1980 1985 1990 1995 2000 2005 2010 2015 
2.0 Year 
1.5 
1.0 
fs 0.5 | m= Before correction 


0.0 | m= After correction 
-1.0 
1975 1980 1985 1990 1995 2000 2005 20102015 
Year 


Fig. 6 Annual RB of precipitation (a), temperature (b), PET (c), maximum temperature (Tmax, d), and 
minimum temperature (Tmin, e) before and after correction in Central Asia during 1975-2014 
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Fig. 7 Temporal variation of water-thermal product index (k index) at annual (a) and seasonal (b—e) scales 
during different periods. The boxes represent the range from the lower quantile (Q25) to the upper quantile (Q75). 
The black horizontal lines represent the means. The upper and lower whiskers extent to the maximum and 
minimum value within the 1.5 interquartile range of the upper and lower quartile, respectively. SSP126, 
SSP1-RCP2.6; SSP245, SSP2-RCP4.5; SSP460, SSP4-RCP6.0; SSP585, SSP5-RCP8.5. SSP, Shared 
Socioeconomic Pathway; RCP, Representative Concentration Pathway. 


3.4 Projected spatial distribution of k index at annual and seasonal scales 


In order to project spatial distribution of k index at annual and seasonal scales under different 
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SSP-RCP scenarios from 2015 to 2100, we selected SSP126 scenario (the minimum emission level) 
as the baseline scenario and used the distribution of k index under SSP126 scenario as the spatial 
benchmark (Fig. 8). Moreover, the difference value of k index between other scenarios (SSP245, 
SSP460, and SSP585) and the baseline scenario (SSP126) was calculated as shown in Figure 9. 

According to Figure 9, we can distinguish the optimal and ordinary scenarios. The optimal 
scenario indicated k index under other scenarios was greater than that under the baseline scenario, 
while the ordinary scenario showed k index under other scenarios was less than that under the 
baseline scenario. In the baseline scenario, the projected k index at annual scale ranged from 
—263.81°C-mm to 43.81°C-mm, showing an explicit decreasing trend from southeast to northwest. 
Higher precipitation, lower temperature, and limited evaporation engendered the relatively high k 
index, accounting for 4.6% of the total area, mostly in the region of Pamir and Tianshan 
mountains; while the negative k index was located in the southern deserts. 

The projected spatial distribution of k index was determined by different seasons and scenarios. 
In terms of seasonality, k index in winter tended to decline from southwest to northeast, which 
was opposite to the spatial pattern of other seasons. Moreover, the variation under SSP-RCP 
scenarios indicated a slightly increased spatial standard deviation of k index in spring, summer, 
and autumn. Furthermore, hydrothermal condition also represented seasonal difference in the 
same region. 

Figures 8 and 9 revealed that, from the perspective of different scenarios, SSP126 scenario was 
the optimal scenario in Central Asia; SSP245 was regarded as the optimal scenario in the western 
Caspian Depression; and SSP460 and SSP585 were the optimal scenario in the eastern plateaus 
and Altay Mountains. From the perspective of seasonality, SSP245 was the optimal scenario in 
the rainfed farming areas of northern Kazakhstan in summer; SSP460 was the optimal scenario in 
the Caspian Depression in autumn; and SSP585 was the optimal scenario in the eastern plateau 
region, Caspian Depression, and Ustyurt Plateau in winter. 


(a) Annual (b) Spring (c) Summer 


À 


(e) Winter (d) Autumn Lengend 


k index (°C-mm) 


geen 
Low: -263181 9__1000km 


ue 


Fig. 8 Projected spatial distribution of k index at annual (a) and seasonal (b—e) scales under SSP126 scenario in 
Central Asia during 2015-2100. Note that this map is based on the standard map (No. GS (2016) 1666) of the 
Map Service System (http://bzdt.ch.mnr.gov.cn/) marked by the Ministry of Natural Resources of the People's 
Republic of China, and the base map has not been modified. 


3.5 Sensitivity analysis of k index to climate variables 


Figure 10 showed the sensitivity analysis results of k index to the four climate variables under 
different SSP-RCP scenarios. The dominant climate variable and its sensitivity were different in 
temporal scales and scenarios. The results indicated that Tmin had little effect on the variation of 
k index (except in winter), whereas k index was most sensitive to Tmax in summer and autumn. 
Spatially, the change of k index was generally sensitive to precipitation in the northern regions, 
while it was sensitive to Tmax in the southern regions. The sensitive ratio of precipitation 
increased from 36.24% under SSP126 scenario to 63.6% under SSP585 scenario, and the ratio of 
Tmax decreased from 57.82% under SSP126 scenario to 28.19% under SSP585 scenario. 

The sensitivity analysis at seasonal scales suggested that the increasing temperatures and 
decreasing precipitation were the main reasons for the declining hydrothermal condition. Under 
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Fig. 9 Projected spatial distribution of the difference value of k index at annual (al—a3) and seasonal (b1—e3) 
scales under SSP245, SSP460, and SSP585 scenarios in Central Asia during 2015-2100. (b1—b3), spring; (c1—c3), 
summer; (d1—d3), autumn; (e1—e3), winter. Note that this map is based on the standard map (No. GS (2016) 1666) 


of the Map Service System (http://bzdt.ch.mnr.gov.cn/) marked by the Ministry of Natural Resources of the 
People's Republic of China, and the base map has not been modified. 


SSP126 scenario, regardless of the spatial difference, the declining hydrothermal condition led to 
a more pronounced hydrothermal contradiction in Central Asia. The projected increase in Tmin 
ranged from 0.04°C-mm/10a to 0.30°C-mm/10a, which was more significant than Tmax and 
temperature. By contrast, the increase in temperature and precipitation was more notable in the 
scenarios with high emission level. In addition, the underlying mechanism of the change of 
hydrothermal condition in the eastern mountains was more complex than that in the rest of 
Central Asia. 


4 Discussion 


4.1 Application and prospect of hydrothermal condition and k index 


Drought is often chosen as a breakpoint in analyzing mismatched hydrothermal condition. Based 
on drought indices, Wang et al. (2021) discovered that severe drought is projected in the southern 
deserts of Central Asia, but its severity is lower than that in the region of Pamir and Tianshan 
mountains. Deng et al. (2020) also found a similar pattern for drought prediction in Central Asia. 
Severe droughts mainly occur in imbalanced regions, whereas regions with relatively balanced 
conditions are more vulnerable to the adverse effects of drought. The spatial distribution between 
the drought degree and k index is similar due to the inherent consistency in hydrothermal 
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Fig. 10 Spatial distribution of dominant climate variable at annual (al—a4) and seasonal (bl—e4) scales under 
SSP126, SSP245, SSP460, and SSP585 scenarios in Central Asia during 2015-2100. (b1—b4), spring; (cl—c4), 
summer; (d1—d4), autumn; (e1—e4), winter. Note that this map is based on the standard map (No. GS (2016) 1666) 
of the Map Service System (http://bzdt.ch.mnr.gov.cn/) marked by the Ministry of Natural Resources of the 
People's Republic of China, and the base map has not been modified. 


condition. Drought indices focus on the disaster strength, but k index is more concerned with 
regional hydrothermal contradictions. Further research is needed to compare the ability of the two 
indices to detect hydrothermal condition in different regions. 

Hydrothermal condition has widespread impacts on vegetation changes at a large scale. 
Previous studies showed that stronger drought triggers vegetation degradation in many regions 
(Horion et al., 2016; Xu et al., 2016; Dubovyk et al., 2019; Yuan et al., 2022). Based on 
multiyear observation data, a study conducted by Ni and Zhang (1997) showed a positive 
correlation between the net primary production of vegetation and k index from taiga to rain 
forest and from temperate desert to deciduous forest. For drought, k index is a valid tool to 
assess the comprehensive hydrothermal states, including balanced and imbalanced conditions. In 
agriculture and livestock industry, k index may have a threshold due to the optimum matching 
condition, the optimum thermal properties under available water resources, or the optimum 
moisture under excessive thermal resources. Therefore, according to the threshold, we can divide 
the k index into different levels, and the projection of hydrothermal condition will be improved. 
Li et al. (2018) attempted to construct a mathematical model between k index and net primary 
productivity, but did not get a reliable result. Further, in areas with higher complexity, such as 
mountains and plateaus, hydrothermal condition is influenced by elevation and terrain (Geng et 
al., 2017; Wang et al., 2022). Future studies need to quantify the influences of terrain by refining 
the spatial resolution in the mountain regions to meet requirements for three-dimensional spatial 
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analysis (Ahmed et al., 2013; Ge et al., 2021). 
4.2 Relationships between hydrothermal condition and agriculture 


We can divide Central Asia into regions with relatively balanced and unbalanced hydrothermal 
condition based on k index, corresponding to rainfed and irrigation agriculture (Kienzler et al., 
2012; Jiang et al., 2019). For example, irrigation agriculture in Uzbekistan and Turkmenistan is 
dominated by unbalanced hydrothermal condition. In contrast, rainfed agriculture and livestock 
industry in the northern Kazakhstan are the result of balanced hydrothermal condition. It is 
suggested that the balanced hydrothermal condition has a positive effect on agricultural 
development. As a part of climate, hydrothermal condition is a limiting but decisive variable for 
agriculture, indicating only the resistance from the natural environment. Some conscious 
interventions, such as irrigation, mitigate the negative impacts of the unbalanced hydrothermal 
condition on agriculture (Zhang et al., 2015; Fang et al., 2019). Furthermore, improved 
agricultural technology and hydraulic engineering reduce resource restrictions, but the impact 
degree is different for water and thermal resources (Reshmidevi et al., 2009; Bannayan et al., 
2010; Carli et al., 2014). An increase in thermal resources in the future period may create a 
development opportunity for agriculture in the Syr Darya River and Amu Darya Rvier. 


4.3 Projected optimal scenario in Central Asia 


According to the results of this study, on the longtime scale, hydrothermal condition is decreased 
under SSP585. Moreover, the unsustainable high emission scenario with a regional superiority 
exerts more hydrothermal pressure on the vulnerable ecosystems and agriculture (Abd-Elmabod 
et al., 2020; Li et al., 2020). In addition, sensitivity analysis suggests that overestimated 
temperatures may lead to a severe hydrothermal contradiction in Central Asia. The 
underestimated precipitation has a similar effect on hydrothermal condition under a low emission 
scenario. The projection of GCMs suggests that reducing the uncertainty of the model is still a 
primary concern for researchers (Guo et al., 2021). 


5 Conclusions 


In this study, we projected hydrothermal condition and analyzed the spatial-temporal variations of 
k index at annual and seasonal scales in Central Asia under different SSP-RCP scenarios during 
2015-2100. In terms of temporal variation, compared with the historical period, projected 
hydrothermal condition is improved from 2046 to 2075 and decreased from 2015 to 2044 and 
from 2076 to 2100. In terms of spatial distribution, the hydrothermal condition is matched in the 
region of Pamir and Tianshan mountains as well as the northern plains of Kazakhstan, and it is 
mismatched in the central and southern deserts. Moreover, the projected hydrothermal condition 
is superior in spring and winter than in autumn and summer. The sensitivity analysis suggests that 
hydrothermal condition in the northern regions is generally affected by precipitation, while 
hydrothermal condition in the southern regions is influenced by Tmax. Further research is needed 
to explore the specific driving mechanisms and the impact of changes in hydrothermal condition. 
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